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ABSTRACT 

Aims. The earliest stages of intermediate- and high-mass star formation remain poorly understood. To gain deeper insights, we study 
a previously discovered protostellar source that is deeply embedded and drives an energetic molecular outflow. 
Methods. The source, UYSO 1, located close to IRAS 07029-1215 at a distance of about 1 kpc, was observed in the (sub)millimeter 
and centimeter wavelength ranges, as well as at near-, mid-, and far-infrared wavelengths. 

Results. The multi-wavelength observations resulted in the detection of a double intermediate-mass protostar at the location of 
UYSO 1. In addition to the associated molecular outflow, with a projected size of 0.25 pc, two intersecting near-infrared jets with 
projected sizes of 0.4 pc and 0.2 pc were found. However, no infrared counterparts to the driving sources could be detected in 
sensitive near- to far-infrared observations (including Spitzer). In interferometric millimeter observations, UYSO 1 was resolved into 
two continuum sources with high column densities (> 10^"^ cm~^) and gas masses of 3.5 M© and 1.2 M©, with a linear separation 
of 4200 AU. We report the discovery of a H2O maser towards one of the two sources. Within an appropriate multi- wavelength 
coverage, the total luminosity is roughly estimated to be ^50 L©, shared by the two components, one of which is driving the molecular 
outflow that has a dynamical timescale of less than a few thousand years. The jets of the two individual components are not aligned. 
Submillimeter observations show that the region lacks the typical hot-core chemistry. 

Conclusions. We find two protostellar objects, whose associated circumstellar and parent core masses are high enough to suggest 
that they may evolve into intermediate-mass stars. This is corroborated by their association with a very massive and energetic CO 
outflow, suggesting high protostellar accretion rates. The short dynamical timescale of the outflow, the pristine chemical composition 
of the cloud core and absence of hot core tracers, the absence of detectable radio continuum emission, and the very low protostellar 
luminosity argue for an extremely early evolutionary stage. 
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1. Introduction 

In the investigation of the earliest stages of intermediate- and 
high-mass star formation, the relative importance of formation 
processes as different as disk accretion and coalescence is still an 
open question (for a recent review, see e.g. B euther et al. 2007 ). 
Detailed multi- wavelength studies are needed to illuminate the 
often confusing observational picture. 

Since high-mass protostars are less ubiquitous than their 
low-mass counterparts, they are on average more distant. The 
resulting low linear resolution often leads to an oversimplified 
picture. This is well illustrated by the study of iBeuther et al.l 
^02a), which reported a case where a single bipolar molec- 
ular outflow discovered previously with single-dish radio tele- 
scopes in a massive-star-forming (MSF) region was resolved 
into a multiple outflow system in interferometric data (see also 
IBeuther etal]l2Q03h . iDavis et al.1 {2004) analyzed near-infrared 
data of two collimated jets in a similar MSF region. They con- 
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elude that these jets are very similar to their low-mass counter- 
parts. 

We report the results of a comprehensive multi- wavelength 
follow-up study of a candidate massive protostellar source that 
was discovered close to IRAS 07029-1215 during a millime- 
ter continuum survey of the surroundings of l uminous IRAS 
sources in the outer galaxy (iKlein et al r i2005h . UYSO 1 is a 
deeply embedded, very young source at a distanc e of only 1 kpc, 
powerin g a h igh-velocity bipolar CO outflow; Forbrich et al.l 
(2Q04b (IPaper D derived the mass of the UYSO 1 cloud core 



from the CO(3-2) line emission to be 44 M©, inferring the hy- 
drogen column density from the integrated CO emission, while 
the mass derived from optically thin dust continuum emission 
(A = 850yum), tracing high column densities, is 15 M©. The 
mass of the outflow was estimated to be 5.4 M©. For details on 
these earlier mass e stimates as well as for a discussion of the 
distance, we refer to "Paper T. Compared to other massive molec- 
ula r outflows, thi s mass is at the lower end of the range reported 
by IBeuther et al.l (l2QQ2bl) . but most of their sources are much 
more distant, at an average distance of 4.7 ±3.4 kpc. No plau- 
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Table 1. Table of relevant positions^ (Epoch 2000) 
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Fig. 1. H2 S(l) emission, as observed with OmegaPrime at the 
Calar Alto 3.5m telescope (filter NB2122), in logarithmic scal- 
ing. The box denotes the close-up VLT view in Fig. O The B 
star HD 53623 is marked by an arrow and the location of IRAS 
07029-1215 is marked by its uncertainty ellipse. 



sible driving source could be identified in IRAS a nd MSX mid- 
and fa r-infrared data. In single-dish observations jB euther et al.l 
(l2008h detected C2H submillimeter emission towards UYSO 1, a 
transition that may be a tracer of the earliest stages of (massive) 
star formation. 

The aim of our follow-up observations of UYSO 1 was to 
study the region at higher angular resolution and to search for 
the driving source of the enormous CO outflow. For these pur- 
poses, we conducted observations in the infrared regime as well 
as in the millimetric to centimetric wavelength ranges. We note 
that for very early evolutionary stages, outflow mass entrainment 
rates and column densities are more conclusive in determining 
whether a massive star forms than a mass determination since 
material is still rapidly accreted. In Section [2l we describe the 
variety of observations carried out before presenting the results 
in Section [3l We discuss our findings in Section (Hand conclude 
in Section [5] 



2. Observations and data analysis 

Near- and mid-infrared observations In the near infrared, ob- 
servations were first performed with the 3.5m telescope on Calar 
Alto, Spain, using the OmegaPrime camera on March 5, 2004, 
both in K' broadband and H2 S(l) narrowband (NB2122). After 
a tentative discovery of two jets, intersecting at the position of 
UYSO 1 , follow-up observations were conducted on several oc- 
casions between January 10 and March 1, 2005 with the Infrared 
Spectrometer and Array Camera (ISAAC) at the Very Large 
Telescope (VLT) of the European Southern Observatory (ESO). 
These observations included Kg, NB2.09 and NB2.13 as well as 
L-band imaging. In January and February 2006, UYSO 1 was 
observed with the mid-infrared instrument VISIR at the VLT 
in the RAH2 filter centered on A = 11.3yum, with an estimated 



Name 



RA (h:m:s) 



DEC (°:':") 



UYSO la 
UYSO lb 
H2O maser 
1st 70 //m peak 
2nd 70 yL/m peak 
24 yt/m peak 



07:05:10.940(1) 

07:05:10.811(3) 

07 05 10.8105(1) 

07:05:10.96 

07:05:11.09 

07:05:11.56 



IRAS 07029-1215 07:05:16.9 



-12:19:00.64(3) 

-12:18:56.84(9) 

-12:18:56.807(3) 

-12:19:09.9 

-12:19:33.1 

-12:19:19.4 

-12:20:02 



^ Numbers in brackets denote the uncertainty in the last digit. 



3cr point-source sensitivity of 9 mJy (using the VISIR Exposure 
Time Estimator v3.2.1). The Calar Alto data were reduced using 
IRAF while we used the ESO Eclipse 5.0 software for the VLT 
data. 



Far-infrared observations We used the Multib and Imaging 
Photometer for Spitzer (MIPS, iRieke et al.ll2004l) onboard the 
Spitzer Space Telescope to obtain maps at 24 fim and 70 jum as 
well as a low-resolution spectrum from 52yum to 105 yum in the 
SED mode. The observations were executed on November, 1 1 
and 12, 2005. The imaging data presented here are post-BCD 
maps processed with version SI 6.1 of the pipeline. Inspection 
of the post-BCD data products and manual reduction with the 
mopex software package showed that the pipeline reduction of 
the data is reliable in our case. The widespread extended emis- 
sion made estimating the background difficult. Some striping 
along the detector columns remain in the 70yum maps. For dis- 
play, the stripes were reduced by smoothing with a three-pixel 
median filter aligned perpendicularly to them. Measurements, 
however, were done on the unsmoothed maps. The MIPS-SED 
data were reprocessed with mopex, starting with the BCD prod- 
ucts created by the pipeline (SI 6. 1.1) to create the spectra along 
the slit. 



Millimeter and submillimeter observations Since the discovery 
observation of the CO outflow, two considerably larger and fully 
sampled CO(3-2) maps of the region were taken on June 4 and 6, 
2005, with the facility receiver B3 at the James-Clerk-Maxwell 
Telescope (JCMT). In December 2003, a map of UYSO 1 and 
its surroundings in C0(2-l) was taken with the Heterodyne 
Receiver Array HERA at the IRAM 30m telescope. However, 
since the CO(3-2) data cover a larger area at a better signal-to- 
noise ratio, we focus our analysis on these data. In July 2005, 
N2H+(l-0), N2D+(l-0), HCO+(3-2) and DC0+(2-l) were ob- 
served towards UYSO 1 with the IRAM 30m telescope and its 
facility receivers in frequency switching mode. 

After first, very short observations with the IRAM Plateau 
de Bure Interferometer (PdBI) in its low-angular resolution D 
configuration in September/October 2003, UYSO 1 was reob- 
served in the newly extended high-resolution A configuration in 
February 2006. While the synthesized beam size in the former 
observations was 13V9 x 5V9 for the 3 mm continuum, this im- 
proved to 2.2'' X 0.6'' in the latter (0.85" x 0.23" at 1 mm). 
Besides the continuum, the CS(2-1) and C0(2-l) transitions 
were observed. In order to be able to account for short-spacings 
in the uv plane, the source was previously mapped in both tran- 
sitions with the IRAM 30m telescope. 

In November 2005, UYSO 1 was observed with the Caltech 
Submillimeter Observatory (CSO) on four occasions (Nov 
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Fig. 2. H2 S(l) emission (NB2.13) in logarithmic scaling, overlaid with the JCMT CO(3-2) observations - total line 10% to 90% in 
10% steps, line wings 30%, 50%, 70% and 90% - as well as the two continuum sources seen with the PdBI (green ellipses indicating 
the FWHM synthesized beam). Four arrows indicate the locations of the jets, and an additional small arrow indicates the eastern 
bow shock of the smaller jet. 



17/18, Nov 19, Nov 20, Nov 22), looking for several submillime- 
ter transitions in the 345 GHz and 220 GHz ranges. The beam 
size£](FWHM) of the CSO are ~ 22'' and -31'' in the 345 GHz 
and 220 GHz bands, respectively. In the same month, the source 
was also observed with the Atacani a Pathfinder Experiment 
(APEX) telescope (Guste n et aPbOOGh to obtain a better C0(3- 
2) spectrum at the source position. The FWHM beam size0 of the 
APEX telescope at 345 GHz is ~ 18". All data from the 30m, 
APEX, and the PdBI were reduced using the GILDAS software 
developed by IRAM and Observatoire de Grenoble. 



^ http : //www. submm.caltech.edu/cso/receivers/beams . html 
^ http://www.apex-telescope.org/telescope/ 



Centimeter radio observations In November 2003, H2O maser 
emission towards UYSO 1 was discovered with the Eff'elsberg 
100m radio telescope. The discovery was confirmed with the 
same telescope in March 2004. At the same time, as well as in 
January 2004, the source was searched for NH3(1,1) emission 
with the Eff'elsberg 100m telescope. On April 13, 2005, it was 
searched for signs of CH3OH maser emission, using the Torun 
32m radio telescope. 

In March 2006, high-resolution radio observations were car- 
ried out with the NRAO Very Large Array (VLA) in A configu- 
ration. Both the 8.4 GHz radio continuum emission as well as 
the newly discovered H2O maser emission were studied. The 
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Fig. 3. Left: continuum- subtracted H2 S(l) emission (NB2.13-NB2.09) in logarithmic scaling. The subtraction of the brightest stars 
left a few artefacts. The box indicates the field of view of the right panel; ellipses indicate the positions of the two millimeter 
continuum sources as in Fig. [2l Right: NIR-continuum-only (NB2.09) close-up view of the central area in linear scaling with the 
3 mm PdBI continuum data shown as contours in steps of 1 mJy. 



Effelsberg data were reduced using the IRAM GILDAS soft- 
ware, and for the VLA data, we used the NRAO AIPS software. 



3. Results 

Based on the new multi- wavelength data, we report the discov- 
ery of two collimated jets that appear to be connected to two pre- 
viously unresolved protostars at the location of UYSO 1. In spite 
of deep searches, we did not detect any infrared counterparts of 
these protostars. Towards one of the two protostars, water maser 
emission was discovered. Observations of molecular transitions 
reveal that UYSO 1 does not show typical hot-core chemistry. 
In the following, we present these results in more detail, diff'er- 
entiating large-scale structure (the clump in which UYSO 1 is 
embedded, the outflow, and the jets) and small-scale structure 
(the protostars UYSO la and UYSO lb) before discussing the 
spectral energy distribution (SED). Relevant positions are listed 
in Tableffl 



3.1. Large-scale structure 

Near- and mid-infrared observations The initial discovery of 
two collimated jets, obtained with the Calar Alto 3.5m telescope 
(Fig.[T]), was clearly confirmed by the VLT observations (Fig. [2]). 
The larger field of view of the Calar Alto narrowband image 
additionally shows H2 S(l) emission in the outer parts of the 
HII region Sh2-297 which is powered by HD 53623. The two 
jets intersect with an angle of 75° at about the submillimeter 
continuum position of UYSO 1, a first indication that UYSO 1 
may harbor more than one source. The larger north- south jet 
has an apparent size of 1.3', corresponding to 0.4 pc at a dis- 
tance of 1 kpc. It also has multiple bow shocks, best seen in the 
continuum- subtracted image in Fig. [3l This jet deviates consid- 
erably from a straight line in its northern lobe, possibly a sign of 



interaction with the surrounding medium. The smaller jet, termi- 
nating in a bright, very compact bow shock at its eastern end, is 
at least 0.2 pc long (it might be even longer, as a barely visible 
chain of features extends well beyond the bright section of its 
western lobe). Both jets have a knotty structure. A quantitative 
estimate of the inclination angles of the two jets with respect to 
the line of sight is difficult. However, the fact that we clearly see 
both lobes of both jets in spite of the surrounding material sug- 
gests that they lie close to the plane of the sky. The narrowband 
luminosities of the two jets, as determined from the 2.12 fim 
data, are L2.i2//m = 0.02Lo for the larger, north-south jet and 
^2.i2//m = O.OOOTL© for the smaller, east- west jet. The luminos- 
ity of the larger flow is at the upper end of the range of values 
found in a survey of Orion A jStanke et al .112002b . a first sugges- 
tion that the driving source is at least of intermediate mass. 

Another prominent feature is the ridge of strong H2 emission 
east of the core, extending from the north to the south, also visi- 
ble in the Calar Alto image. This ridge appears to be the edge of 
the molecular clump where UYSO 1 is embedded in since also 
the CO and millimeter continuum emission drop dramatically in 
that region. The emission ridge indicates where the H II region is 
interfacing the molecular cloud. Whether the H2 emission is ex- 
cited by shock fronts in the ionisation front or by UV pumping in 
the photon-dominated region (PDR) between the H II region and 
the molecular cloud, we cannot decide, although a PDR appears 
more likely (see the discussion of oxygen fine structure lines in 
Section [33]) . This would require th e comparison of the relativ e 
intensities of several H2 lines (e.g. jHoUenbach & Shull|[T977h . 
data that are not available at present. The H2 emission in jets is 
usually shock-excited. 

Far-infrared observations Neither 24yum nor 70yum mid- 
infrared sources were found to be in direct relation to the mil- 
limetre continuum peak (Fig.©, even though the 70yum obser- 
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Fig. 4. MIPS 24yum (left) and 70yum (right) data. Upper panel: original maps, in logarithmic scaling. Contours of the data (grey) 
and of the CO(3-2) emission (black) are displayed as reference. Lower panel: maps with extended emission filtered out (see text), 
in linear scaling. Contours of the CO(3-2) emission are displayed as reference. The color scale is diff'erent from the one used in the 
upper panel. In all plots, the locations of UYSO la and UYSO lb are marked with small ellipses, as in Fig.[2l 



vations may show some emission at that position (see discussion 
of the SED in Section [33l) . Instead, the 24yum emission closely 
follows the H2 S(l) emission ridge as the PDR is also heating the 
dust at the surface of the clump. The TOyum emission also starts 
at this rim, but extends further into the cloud as it traces colder 



material deeper in the cloud. The emission peak in the 70 fim im- 
age is separated by 13'' to the north-west from the 24yum peak. 
It is still another 11'' away from the millimeter sources. In Sect. 
3.3, we discuss the possibility of UYSO 1 being swamped by 
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Fig. 5. CO(3-2) spectrum towards UYSO 1, as observed with the 
APEX telescope, in [K]. In order to show the line wings, only 
the lower 10% part of the line are shown. The grey line shows a 
Gaussian profile fitted to the entire line. 
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Fig. 6. The two 3 mm continuum observations using the PdBI, at 
the same angular scale (grey scale and solid contours: A config- 
uration, contour step: 1 mJy; dashed contours: D configuration, 
contour step: 2 mJy). The brighter source is UYSO la, the fainter 
one is UYSO lb. 



surrounding extended emission and try to determine its potential 
contribution to the FIR emission. 

(Sub-)millimeter and centimeter radio observations The new 

CO(3-2) single-di sh map, shown in Fig. [2l corroborates the re- 
sults presented in iPaper U There is a strong gradient in emis- 
sion towards the H II region east of the position of UYSO 1 . The 
APEX CO(3-2) spectrum taken at the position of UYSO 1 shows 
the line wings most clearly (Fig.O. The outflow is prominent in 
both the CO(3-2) and the C0(2-l) maps. Within the uncertain- 
ties, the molecular outflow coincides with the larger of the two 
NIR jets. Thus, in the sin gle-dis h CO maps, one massive out- 
flow already discussed in Paper t is detected, but the new, larger 
CO(3-2) map additionally shows an indication of an east- west 
outflow, at least at its western end. Possibly, the millimeter out- 
flow related to the smaller jet is weaker than the large previously 
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Fig. 7. H2O maser emission towards UYSO 1, as observed with 
the Efl'elsberg 100m telescope, axis units are Jy. The VLA de- 
tection of an H2O maser in UYSO lb has an unresolved velocity 
of 23.3 km s"^ marked by an arrow, with a peak flux density of 
0.92 Jy. The system velocity is v=12 km s"^ . 



known outflow and as such difficult to detect when both systems 
are superimposed. 

From the CO(3-2) map, the mass of the outflow can be esti- 
mated in direct comparison to Paper I, following the procedure 
outlined there as well as in Henning et al.l (|2000), using the pro- 
portionality of the integrated CO main-beam temperature and 
H2 column density. The new, larger map has a slightly smaller 
SNR than the data used in Paper I, afl'ecting the line wings. For 
the entire line profile, the 40% intensity contour (which is not 
entirely within the map) traces 33 M©, the red- and blueshifted 
line wings trace 2 M© and 1 M©, respectively (down to the 10% 
contour). These numbe rs are sli ghtly lower than those derived 
from the deeper data in iPaper 1 where the cloud mass was es- 
timated to be 40 M©. While it is difficult to give quantitative 
uncertainties for these mass estimates, we note that the masses 
are uncertain by at least a factor of a few, i.e., they are order- 
of-magnitude estimates. While for the extent and the maximum 
velocities of the outflows, the original JCMT data in the outflow 
lobes are still the best due to their SNR, we note that our deeper 
APEX pointing on UYSO 1 shows line wings that are compara- 
ble to what was previously only detectable in the ou tflow lobes. 
The outflow velocity Vproj = 30 km s"^ from lPaperl compares to 
Vproj = 26.5 kms"^ derived at the position of UYSO 1 (Fig. [51 
compare to Fig. 4 in Paper I). In the new, large CO(3-2) map, the 
10% contours of the line wings indicate a total size of the molec- 
ular outflow of 52", or 0.25 pc. Compared to Paper I, the inferred 
size of the outflow is basically unchanged. Assuming that we 
have an edge-on disk geometry rather than a pole-on view (as is 
suggested by the clear detection of both near-infrared jet lobes) 
allows us to better constrain the outflow and its mass entrainment 
rate. An inclination angle of 80° with respect to the line of sight 
yields an outflow dynamical timescale (td(i) = /?out(0/vout(0 of 
only a few hundred years and an outflow mass entrainment rate 
of M = 4 X 10"^ M©yr"^. Even for an inclination of / = 57.3° 
(see Paper I), the dynamical timescale is still only 3000 years 
with a mass entrainment rate of M = 1 x 10"^M©yr"^ We 
can conservatively constrain the dynamical timescale to less than 
10"^ yr except for very small inclination angles (< 25°) and the 
mass entrainment rate correspondingly to > 3 x 10"'^M©yr"^: 
as discussed above, also the NIR observations of the jet lobes 
suggest large inclination angles with respect to the line of sight. 
Also, in a pole-on configuration with small inclination angles, 
we should see the driving sources in the NIR/MIR when looking 
into the outflow cavities, which is not the case. Based on these 
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Table 2. Molecular transition lines observed towards UYSO 1 



Transition 


Telescope 


T*a [K] 


iMrl3(^i,i ) 


rLrreisDerg luum 


< U. i J Jv 


N2H+(l-0) 


IRAM 30m 


< 0.1 K 


HCO+(3-2) 


IRAM 30m 


3.5 K 


DC0+(2-l) 


IRAM 30m 


0.5 K 


''CO 


CSO^ 


1.8 K 


HCN(4-3) 


CSO 


0.7 K 




CSO 


< 0.4 K 


H2CO(5(i,5)-4(i,4)) 


CSO 


0.3 K 



" notable undetected molecules in our 345 GHz-range CSO 
observations: CH3CN, CH3OH, SO, SO2, HNCO, HCOOCH3, 
HCCCN, CH3CH2CN, and additionally in the 220 GHz range: SiO 
^ blends at 349.338 GHz and 349.400 GHz 



outflow properties, we estimate and discuss the accretion lumi- 
nosity in Section m 

New molecular line observations beyond those in lPaper ll are 
summarized in Table [2l To better trace the o verall mass, ^ ye als o 
studied UYSO 1 in 1^00(2-1). Following Scovi Ue et alJ (Il986h . 
the emission peak at the position of UYSO 1 corresponds to 
a mass of 8 M© in the CSO beam size of ~ 33" (FWHM). 
Interestingly, neither NH3 nor N2H^(l-0) emission was found 
towards the source, with upper limits of <0.15 K and < 0.1 K, 
respectively. HCO+(3-2) and DC0+(2-l) were detected. The 
HCO^(3-2) line may consist of several velocity components. 
Several submillimeter molecular lines that we searched for using 
the CSO were not detected (see footnote in Table O. These re- 
sults indicate that UYSO 1, with only very few detectable molec- 
ular transition lines, does not show typical hot-core chemistry. 

Based on the upper limit of N2H^(l-0) emission, we can es- 
timate corresponding upper limits for the column density and the 
abundance. F or the column density, A^a^2^+ ~ 8 x IO^^AvTr cm"^ 
(iBenson et a l. 1998), we derive an upper limit of Nn^h+ ~ 1.6 x 
10^^ cm"^, assuming a line width of 1 km s"^ . An upper limit for 
the abundance follows when relating this to the hydrogen col- 
umn density derived by lKlein et al.l (l2005b from the submillime- 
ter radio data at 850 yum, Nh = 5.8 x 10^^ cm'^ (for T = 20 K). 
The result, an abundance of only Nn^h+/Nh < 2.6 x 10"^^, is 
surprisingly low compared to values of ^ 10"^^ found in starless 
cores bv .Tafalla et al.l (l2004 . 

3.2. Small-scale structure 

Millimeter observations The millimeter continuum data from 
the two PdBI observations are ideal to study the small-scale 
dust continuum emission. In the low-resolution (D configura- 
tion) data, only a single source is detected at 3 mm. In the high- 
resolution data (extended A configuration), this source is re- 
solved into two protostars, UYSO la and UYSO lb (Table [B; 
both sources remain unresolved. At 1 mm, only the brighter 
southern source, UYSO la, is clearly detected due to limited sen- 
sitivity. Results of the two 3 mm continuum datasets are shown 
inFig.[6l 

We follow iBeuther et aP (I2003L l2005h in determining the 
mass from the 3 mm continuum data, as traced by the optically 
thin dust emission of the proto st ellar envelopes. We assume a 
temperature of T = 45 K (iPaper ih while leaving t he gas-to-dust 
ratio (100:1) and other quantities the same as in iBeuther et al.l 
(l2QQ3h who estimate the results to be accurate within a factor of 
five. The single source detected in the low-resolution data cor- 
responds to a gas mass of 9.5 M©. In the high-resolution data. 



two continuum sources are detected in the 3 mm band. These 
are at roughly the positions of two faint near-infrared contin- 
uum sources, but do not fit exactly (see discussion of the NIR 
observations below). Both objects appear to be intermediate- 
mass protostars with gas masses of 3.5 M© for the more mas- 
sive component, UYSO la, and 1.2 M© for the second compo- 
nent, UYSO lb. The linear separation between the two sources 
is 4V17, or 4200 AU. The hydrogen column densities corre- 
sponding to the derived gas masses are A^h = 6.2 x 10^^ cm"^ 
and A^H = 2.2 x 10^4 cm'^ for UYSO la and UYSO lb, re- 
spectively. As a crude approximation, we convert these to vi- 
sual extinction according to the empirical relation A^H[cm"^] ^ 
2 X 10^^ X Av[mag] (Ryter 1996; Vuong et al. 2003), resulting 
in visual extinctions of Ay >1000 mag. We note that the de- 
rived column densiti es are well above typ ical values for low- 
mass protostars (e.g. jMotte & Andril2QQll) and also above the 
minimum column de nsity o f 3 x 10^^ cm"^ (or Igcm"^) that 
Krum holz & McKe3 (l2QQ8l) derived for the formation of mas- 
sive stars. 

The molecular line data collected with the PdBI suff'ers from 
missing flux, and the two extreme configurations are difficult to 
combine due to virtually no overlap in the uv plane. Thus, we 
only briefly discuss the CS(2-1) data here. While in the D con- 
figuration, again only a single source is detected with barely not- 
icable velocity structure, the A configuration data show emission 
that is largely resolved out. 

Centimeter radio observations H2O ma ser emission, a sign - 
post of low- and high-mass star formation (iHenning et al.lll992h , 
was discovered in single-dish observations towards UYSO 1 in 
two velocity components symmetrically spaced around the sys- 
tem velocity of visr=12 kms"^ (Fig.|3- One of the two compo- 
nents has a double-peaked substructure. In high-resolution VLA 
A-array observations, however, only a single maser spot was 
found at the position of the north-western millimeter continuum 
source, UYSO lb, close in velocity to the brighter component 
in the Eff'elsberg maser spectrum. No CH3OH maser emission 
was detected (< 0.3 Jy). In the centimeter continuum at wave- 
lengths of 1.3 cm and 3.5 cm, only 3cr upper limits of 5.6 mJy 
and 0.21 mJy for the flux densities could be determined, respec- 
tively. Thus, no detectable H II region has formed yet. 

Near- and mid-infrared observations In deep NIR imaging, 
two sources were detected in the NB2.09 continuum and the L 
band which are close to the positions of the millimeter proto- 
stars, but they do not coincide with them (Figs. [2l and [3]). The 
astrometric accuracy of the near-infrared data, refined with po- 
sitions from the 2MASS catalogue, is estimated to be < 0.2''. 
Given the enormous visual extinction towards the two proto- 
stars (see above), \ ye prob ably only se e scattered ligh t from 
their vicinities re.g. JLinz et al. 2005; Wei gelt et all2006h . Since 
UYSO la lies in the large north- south jet, and UYSO lb lies in 
the east- west jet, the two millimeter sources probably are their 
driving sources. Then, UYSO la would also power the dominant 
molecular outflow coinciding with its jet while UYSO lb excites 
the water maser. Notably, the NB2.09 NIR continuum source just 
to the south-east of UYSO la and the source just to the east of 
UYSO lb appear not to be point sources; it is also noteworthy 
that both appear on the sides of the blueshifted outflow lobes, 
which would be tilted towards us and thus sufl'er less extinction. 
In the VISIR observations at 1 1 yum, no source was detected at 
the position of UYSO 1. As an upper limit, we use the expected 
3cr point source sensitivity of 9 mJy. 
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3.3. Spectral Energy Distribution 

Compa red to the spectral energy distribution of UYSO 1 in 
iPaper ll which was based on the submillimeter detections and 
upper limits derived from IRAS data, we now have much more 
comprehensive information to constrain the combined luminos- 
ity of UYSO la and lb. 



Broad-band Spitzer observations As noted in Sect. [3711 
UYSO 1 remains undetected in the MIPS scan maps. Bright ex- 
tended emission from the clump around UYSO 1 is detected but 
no compact emission peak that could be ascribed to the core con- 
taining UYSO la and lb is seen. The question is how to estimate 
upper limits for the flux density of UYSO 1 at 24 jim and 70 jim. 
As a conservative estimate, we perform aperture photometry at 
the position of UYSO 1, using apertures of W and 32'' in diam- 
eter with background annuli from 80'' to 100" and 78" to 130", 
for 24 and 70 yum, respectively. These apertures include UYSO 1 
but do not cut into the FDR, whereas the background annuli lie 
beyond the main extended emission. We applied the appropri- 
ate aperture and color corrections as described in the MIPS Data 
Handbook v3.3 as derived from the theoretical PSF. This method 
yields conservative upper limits of 0.3 Jy and 70 Jy for UYSO 1 
at 24 Jim and 70 yum, respectively. 

However, these values mainly measure the large-scale emis- 
sion of the externally heated envelope at the position of UYSO 1 . 
Inspired by unsharp masking, we estimate the extended emis- 
sion by convolving the maps with a Gaussian profile and then 
subtracting it to retain the small-scale structure. The FWHM of 
the Gaussian has been 40". The unsharp-masked maps are dis- 
played in the lower panels of Fig. (H The small-scale structure 
is now much more pronounced, but no peak at the position of 
UYSO 1 becomes apparent. Repeating the aperture-photometry 
should now result in values much less aff'ected by the extended 
emission. Virtually no flux is left in the 24 yum map at the posi- 
tion of UYSO 1 . Only an upper limit of 0. 1 Jy can be derived. At 
70 yum, the aperture photometry retains 25 ± 2 Jy at the position 
of UYSO 1. This number is a second upper limit since there is 
no compact emission apparent at this position (although there is 
a conspicuous extension of the 70 jim emission towards the posi- 
tions of UYSO la and b, possibly indicating that these protostars 
start to become visible at about this wavelength). There is con- 
siderable uncertainty as to how well the extended emission and 
the cloud background are subtracted. 

Low resolution Spitzer spectroscopy (SED mode) Even 
though UYSO 1 remains undetected in the MIPS imaging data, 
the low-resolution spectroscopy can help characterizing the sur- 
roundings. The inset of Fig. [8] shows the location, the size, and 
the pixels of the MIPS slit relative to the 70 yum imaging data; 
the resulting spectrum is shown in the main part of the figure. 
This spectrum was extracted from the three pixels closest to the 
position of UYSO 1. The cloud background was removed by 
subtracting the average of the three pixels east and three pix- 
els west of the UYSO 1 spectrum. All the spectra are derived 
by averaging over three pixels along the slit and applying the 
proper aperture correction. The 3cr error bars are the uncertain- 
ties propagated from those reported in the post-BCD products. 
The spectrum derived for the position of UYSO 1 is compatible 
with the flux derived from the 70yum map. It does not show any 
strong features and the slope can be described by a black body 
of temperature 40 K. Note again that UYSO 1 itself is probably 
not detected. 
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Fig. 8. MIPS Spectrum of UYSO 1 with 3cr error bars; in grey a 
40 K black body spectrum; inset: Location, size, and integrated 
intensity of the MIPS slit on 70yum-map with 850yum-contours 
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Fig. 9. The [Oi] and [Oiii] line: The intensity in the MIPS slit 
are the [Oi] 63yum-line on top of the 70yum-map (850yum con- 
tours). The inset displays again the variation of the [Oi] 63 yum- 
line across the slit plus the [Oiii] 88yum-line (grey). Error bars 
are 3cr. 



Oxygen fine structure lines A feature that is detected in the 
Spitzer SED-mode observations despite its low spectral reso- 
lution are oxygen fine structure lines at the eastern rim of the 
molecular core. Fig.[9]again shows the 70 yum map with the MIPS 
slit, this time with the intensities of the [Oi] line at 63 yum. The 
inset shows the variation of the [Oi] line and the [Oiii] (88yum) 
line along the slit. The subtracted continuum was obtained by 
averaging the two pixels neighboring the line-containing pixel 
in wavelength. 

We clearly see here a PDR radiating in two important cool- 
ing lines on the eastern side of the core where the radia- 
tion from IRAS 07029-1215 impacts the molecular cloud. The 
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emission comes only from a small area as atomic oxygen is 
only present up to optical depths of Ay ~ 10 mag, and the 
PDR must be seen edge on. The [Oi] line's intensity of a few 
10"^ erg cm~ ^s~^sr~^ is expected for PDR s with densities of 10^ 
to lO^cm-^ (iHoUenbach & Tielenslfl999h . 



4. Discussion 

In view of our previous knowledge on UYSO 1 , the main new re- 
sult that warrants discussion is its newly constrained SED. Also, 
the nondetections in N2H^ and NH3 are of special interest. 

The non-detection of UYSO 1 at mid-infrared wavelengths is 
surprising given the earlier estimate of the SED in Paper I. The 
mid-infrared part of that SED was only weakly constrained by 
upper limits from IRAS data, clearly contaminated by the nearby 
bright source IRAS 07029-1215. The previous estimate includes 
the entire cloud core, heated also externally. The previous lumi- 
nosity upper limit of < 1900Lo corresponds to the large-scale 
emission in the Spitzer maps. Between 24 and 70yum, the lumi- 
nosity of the filtered-out large-scale emission is ~ 1500 L©. 

We can estimate the amount of energy that is intercepted and 
re-radiat ed by the PDR neighboring UYSO 1 . The luminosity es- 
timate in lPaper 11 includes any such emission inside the diameter 
of 80% encircled energy of IRAS at 100 jim. Our MIPS data in- 
dicate that extended emission is indeed present around UYSO 1 
on these scales. The projected distance between UYSO 1 and 
HD 53623 is 0.46 pc. Since the cloud is illuminated from the 
side, the real distance should not be much different, i.e., <1 pc 
(if it protrudes less than 60° out of the plane of the sky). When 
seen from HD 53623, the size of the extended emission sub- 
tends an angle of 28-60°, or a solid angle of 0.18-0.85 sr, 
assuming that the cloud is at a distance of 1-0.46 pc and ap- 
pears circular. The remaining uncertainty is the luminosity class 
of HD 53623 for which spectral types of BIV (Claria 1974) 
and BlII/III (Houk & Sm ith-Moordll988l) are given in the lit- 
erature. C orrespondingly assuni ing luminosities of I6OOOL0 or 
39000 Lo (ISchmidt-Kaleilll982h . the cloud core intercepts and 
eventually re-radiates 230-2640 L©. With the uncertainties in- 
volved, it seems reasonable to assume that the previously esti- 
mated upper limit for the luminosity was dominated by exter- 
nally heated extended emission. 

In spite of the new data, the SED is only constrained by up- 
per limits in the infrared regime. Since the newly discovered two 
components are not resolved in the submillimeter continuum, we 
use the unresolved PdBI continuum data at 3 mm (from the D 
configuration) and do not take into account the high-resolution 
PdBI data for the SED fit. Assuming an isotropic radiation field, 
a new modified-blackbody fit (with dust opacities of ATy oc and 
P = 2) yields a luminosity upper limit of ~ 50 L© for both com- 
ponents combined (see Fig. [T0|), assuming a distance of 1 kpc. 
This may be a severe lower bound if most of the luminosity es- 
capes along the outflow cavities and is not re-processed into in- 
frared radiation. 

In spite of the uncertainties involved, it is interesting to com- 
pare this luminosity to an estimate of the accretion luminosity. 
The outflow mass entrainment rate translates into a mass loss 
rate of the driving jet which in turn translates into the actual ac- 
cretion rate onto the star. Beut her et al.l |2002b) estimate that this 
actual accretion rate is lower than the outflow mass entrainment 
rate by a factor of about 6. The above-mentioned conservative 
lower limit for the mass entrainment rate thus translates into a 
continuous accretion rate of > 5 x 10"^ Moyr"^ A mass entrain- 
ment rate of M = 4 x 10"^ M© yr"^ translates into an accretion 
rate of 7 x 10""^ M© yr"^ . Based on the accretion rate, it is possible 



to estimate the accretion luminosity. Lace = G • • Mace • • Of 
course, the stellar parameters are barely constrained in this case, 
and the situation is further complicated by the assumption of 
isotropic radiation, but we can carry out an order-of-magnitude 
chec k. For a 3 M© star and a corresponding radius of 5-8 Rq 
(Pal la & Stahleilll992l: [Yorke & Bodenheimeill2008h . the lower 
limit of the accretion rate corresponds to an accretion luminosity 
of 600-900 L©. For comparison, a 10 M© star with a correspond- 
ing radius of 7-12 Rq would have an accretion luminosity of 
1300-2200 L©. The estimated accretion luminosity thus appears 
to be about an order of magnitude larger than the luminosity up- 
per limit deduced from the SED with the assumption of isotropic 
radiation. However, we note again that, given the likely presence 
of accretion disks, the radiation field is highly anisotropic. In 
particular, such disks would have the highest extinction, poorly 
determinable in the mid-infrared, in the direction towards the 
observer. Therefore, and due to the fact that the accretion lumi- 
nosity is but a simple estimate with uncertain assumptions, the 
above discussion should be regarded as very tentative. 

It remains unclear whether the empirical relation be- 
tween the mass entrainment rate of the outflow andjthe lu- 
minosity of the central object t hat wa s used in iPaperll is 
rehable dShegherd & Churchwelll 119961: iHenning et al.l I200QI: 

etaDE 



iBeuther et al.H2002bl see also IWu et al.ll2005h . We note never- 
theless that the mass entrainment rate appears to be higher than 
what would be expected when applying this relation for a driving 
source with a luminosity of -50 L©. 

UYSO 1 was not observed to show centimetric NH3 emis- 
sion or millimetric N2H ^ emission alth ough we detect HC0^(3- 
2) at the same position. FWomack etal.l (L1991) find a similar ef- 
fect towards Orion-IRc2 and speculate that N2H^ may be de- 
pleted b}^_shocks_o£_wmds_^^ to that source. On the other 
hand jTafalla & Santiago! (l2004l) found an extremely young star- 
less core with constant C^^O and an unusually low N2H^ abun- 
dance. The fact that we do not see typical hot-core chemistry 
towards UYSO 1, combined with the still low luminosity of the 
two sources suggests that we see very young sources. While it 
remains unclear how exactly the outflow fits into this picture due 
to the unknown inclination angle, a sufficiently large inclination 
angle would explain the extinction towards the central source 
and lead to a high mass entrainment rate of the outflow. 

5. Conclusions 

We report the detection of two highly collimated jets intersect- 
ing close to the position of the previously identified candidate 
massive protostar UYSO 1. Follow-up high-resolution millime- 
ter radio observations show two continuum sources with masses 
of 3.5 M© (UYSO la) and 1.2 M© (UYSO lb) and high column 
densities (> 10^"^ cm"^) close to the geometric centers of the jets 
which probably are their powering sources. The projected dis- 
tance between the two sources is 4200 AU. UYSO la appears 
to power the energetic molecular outflow that was previously 
discovered towards UYSOl. UYSO lb contains a previously 
unknown H2O maser source while a search for CH3OH maser 
emission remained unsuccessful. The two millimeter protostars 
are not detected in deep near-infrared imaging and remain unde- 
tected also at mid- to far-infrared wavelengths. Also, no devel- 
oping H II regions were found. Submillimeter observations show 
that the region is rather cold and does not show typical hot-core 
chemistry. Curiously, not even emission in NH3 or N2H^ was 
detected. An attempt to constrain the combined SED of the two 
sources yields an estimated luminosity of roughly L ^ 50 L©, 
which is about an order of magnitude lower than a tentative es- 



10 



Forbrich et al.: A multi- wavelength study of a double intermediate-mass protostar 



10* r- 
10^ r 

10^ r 



10' r 




1 10^ 10^ 10^ 10* 10^ 

Fig. 10. The SED of the combined two sources that make up 
UYSO 1. The fit mainly relies on the submillimeter continuum 
data and the Spitzer upper limits and thus again constitutes only 
an upper limit in its infrared part. The two data points at 3 mm 
correspond to the high-resolution and low-resolution PdBI data. 
Fitting modified-blackbody radiation to this SED yields an upper 
limit for the estimated total luminosity of -50 L©. 



timate of the accretion luminosity. The key to combining these 
new results into a coherent picture probably lies in the unknown 
outflow inclination angle with respect to the line of sight. A 
large inclination angle would lead to configurations in which 
the circumstellar disk is seen nearly edge-on, obscuring the cen- 
tral object and helping to explain their infrared non-detections. 
For moderately large inclination angles of about > 25 °, where a 
thick circumstellar disk and a dense envelope would obscure the 
central object, the molecular outflow powered by UYSO la has 
a dynamical timescale of < 10000 years with enormous outflow 
mass entrainment rates of > M = 3 x lO'^Moyr"^. The lumi- 
nosity of the corresponding near-infrared jet is comparable to the 
most luminous jets found in a survey of Orion A. In addition to 
the peculiar chemistry and the still relatively low luminosity, this 
indicates a very early evolutionary stage. 
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